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Abstract 

Two-dimensional periodic surface waves propagating under the combined influence 
of gravity and surface tension on water of finite depth are considered. Within the 
framework of small-amplitude waves, we find the exact solutions of the nonlinear 
differential equation system which describes the particle motion in the considered 
case, and we describe the possible particle trajectories. The required computations 
involve elliptic integrals of the first kind, the Legendre normal form and a solvable 
Abel differential equation of the second kind. Some graphs of the results are included. 

1 Introduction 

We are concerned with the motion of periodic plane waves which are propagated 
in water of finite depth and with free surface, under a gravitational field. The 
governing equations of motion are the incompressible Euler equations. On the 
free surface we take into account the influence of surface tension. The surface 
tension will appear in the formulation of the boundary conditions. Thus, one 
obtains the boundary- value problem for capillary-gravity waves. There are very 
few explicit solutions known for the water-wave problems. The first such solu- 
tion for pure gravity water waves, was described by Gerstner |19j . This solution 
was independently rediscovered later by Rankine [37] . Modern detailed descrip- 
tions of this wave are given in recent papers [2] and [23]. Gerstner's solution 
is restrictive: it exists only for deep water and it describes a rotational wave. 
Related to Gerstner's solution, Constantin, in [3], constructed an explicit ro- 
tational solution to the nonlinear governing equations describing gravity water 
waves that progress along the shoreline. Beneath Gerstner's waves it is possible 
to have a motion of the fluid where all particles describe circles with a depth- 
dependent radius ([2], [23] ). Crapper [14] derived an exact solution for pure 
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capillary waves travelling at a constant velocity at the surface of a fluid of infi- 
nite depth. He showed that as the amplitude increases, the waves develop broad 
crests and sharp troughs. Crapper's solutions ultimately reach a limiting form 
for which there is a trapped bubble at the troughs. For higher values of the am- 
plitude the free-surface profiles are self-intersecting and therefore non physical. 
In [24 , by the use of the general method for calculating trajectories given by 
Longuet-Higgins [34], the particle trajectories in Crapper's nonlinear capillary 
waves are derived. It is found that the orbits of the steeper waves are neither 
circular nor closed. Crapper's solution was extended in the case of a fluid of 
finite depth by Kinnersley [35] . The particle trajectories for Kinnersley's waves 
were calculated in [26] . For nonlinear capillary-gravity waves no exact analytic 
solution has yet been found. Making use of numerical studies, in [25] the parti- 
cle trajectories in irrotational nonlinear capillary-gravity waves are investigated. 
There are several unanswered questions about water waves with surface tension 
and vorticity. The existence of regular periodic travelling waves with vorticity 
was recently established (see [TT], [H]). For steady periodic gravity waves the 
symmetry is known to be ubiquitous (see [B], [H]). The study of the symme- 
try of rotational water waves was initiated in [7], [8]; for irrotational flows see 
also [36j . However, exact information about the flow beneath such waves is not 
readily available, even in the irrotational case. This paper addresses this issue. 

It is fortunate that a great number of observations can be explained on the 
basis of the small-amplitude wave theory. Within this framework, in this pa- 
per we investigate the capillary-gravity waves and the internal motion of the 
fluid under the passage of such waves. We simplify the full system of equa- 
tions by a linearization found in [28], which is around still water and which is 
slightly different from the classical case in line with the Stokes condition for 
irrotational flows (for the latter, see for example, [5J and [23] where pertur- 
bations of laminar flows are considered, that is, flows characterized by a flat 
surface r\ = 0; these two approaches to linearization are surveyed in [18]). Af- 
ter rewriting the governing equations of motion and the boundary conditions 
for the capillary-gravity water-wave problem in an appropriate non-dimensional 
form, three non-dimensional parameters arise: e =amplitude/depth, the am- 
plitude parameter, S =depth/wavelenght, the shallowness parameter and W e , 
called Weber number, which comes from the surface tension on the free surface. 
These parameters help us to characterize the various types of approximation. 
We suppose that the water flow is irrotational; thus, in addition to the full sys- 
tem of equations we also have the irrotational condition that we write in the 
suitable set of non-dimensional variables. The linearized problem is obtained 
by letting e — > 0, S and W e being fixed. Solving this problem, we obtain a pa- 
rameter Co by which we can describe different backward flows in the irrotational 
case: still water (cq = 0) , favorable uniform current Cq > and adverse uniform 
current cq < 0. After getting the general solution of the linearized water-wave 
problem we investigate further the nonlinear equations of the motion of the fluid 
particles to offer much insight into the fluid motion. 

It was widely believed that as the small-amplitude waves propagate on the 
surface the particles of the fluid move on closed orbits. Analyzing the first-order 
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approximation of the nonlinear ordinary differential equation system which de- 
scribes the particle motion, it was indeed obtained that all water particles trace 
closed, circular or elliptic, orbits (see, for example, [TS], [3T], [33], [35], [3J5], [30]V 
But in |10) it was proved, using phase-plane considerations for the nonlinear or- 
dinary differential equation system which describes the particle motion, that in 
linear irrotational gravity water waves no particle trajectory is actually closed, 
unless the free surface is flat. Each particle trajectory involves over a period a 
backward/forward movement, and the path is an elliptical arc with a forward 
drift; on the flat bed the particle path degenerates to a backward/forward mo- 
tion. In [TH] recent results in the understanding of particle paths within different 
types of progressive water waves are surveyed, in the framework of linear theory 
as well as in the framework of exact theory of periodic symmetric waves, and 
in the presence or not of the background currents and vorticity. The results are 
in agreement with Stokes' observation for a net mass drift jJT]. In the linear 
framework, by using phase-plane considerations, one obtains that the particle 
trajectories in linear deep-water waves (see in linear gravity water waves 
over finite depth with constant vorticity (see [IB], [H]), and in linear irrotational 
capillary and capillary-gravity water waves (see [21], [22]) are not closed. Small- 
amplitude shallow-water waves were studied in [55] (the irrotational case) and 
[2"9"] (the constant vorticity case), and the exact solutions of the nonlinear ordi- 
nary differential equation system which describes the particle motion below such 
waves were found. Depending on the strength of the underlying uniform current 
[28] or depending on the relation between the initial data and the constant vor- 
ticity [29] , it was obtained that some particle trajectories are undulating curves 
to the right, or to the left, others are loops with forward drift, or with backward 
drift, and others can follow some peculiar shapes. In the framework of exact 
theory, the forward drift of the particles as the wave progresses was given by an- 
alyzing a free-boundary problem for harmonic functions in a planar domain (see 
[4] for Stokes waves, [9] for solitary waves and [20] for deep-water Stokes waves) 
or by applying local bifurcation theory (see [J3] for small-amplitude waves with 
vorticity) . 

In this paper we continue the study started in [30 . We provide explicit 
solutions for the ordinary differential equation system describing the motion of 
the particles beneath small-amplitude capillary-gravity waves which propagate 
on the surface of an irrotational water flow with a flat bottom. In the case 
when the constant en equals the non-dimensional speed of propagation of the 
linear wave, the required computations involve elliptic integrals of the first kind 
and their Legendre's normal form. The six exact solutions obtained in this case 
contain Jacobian elliptic functions in their expressions. In the case when the 
constant cq is different from the non-dimensional speed of propagation of the 
linear wave, the computations involve a solvable Abel differential equation of 
the second kind. In both cases we remark that the solutions obtained are not 
closed curves. We analyze some solutions in detail and we draw their graphs. 



2 Preliminaries 



4 



2 Preliminaries 

2.1 The water-wave problem 

For two-dimensional periodic waves the motion is identical in any direction 
parallel to the crest line and is periodic in the spatial direction in which the wave 
is propagating. To describe these waves we consider a cross section of the flow 
that is perpendicular to the crest line with Cartesian coordinates {x,z), the x- 
axis being in the direction of wave propagation and the z-axis pointing vertically 
upwards. The water flow under consideration is bounded by a rigid horizontal 
surface below at z — and a free surface above at z — h + ri(x,t), where the 
constant h > is the mean water level. We denote by (u(x, z, t), v(x, z, t)) the 
velocity of the water. Let gravity act now in concert with surface tension. The 
surface tension will play a role in the formulation of the boundary conditions 
but not in the equations of motion valid in the fluid domain. Assuming that 
the water is both homogeneous (constant density p) (see [35]) and inviscid, we 
obtain within the fluid domain the equation of mass conservation (MC) together 
with Euler's equations (EEs) (see [31]). The boundary conditions for the water- 
wave problem are the kinematic boundary conditions as well as the dynamic 
boundary condition. The kinematic boundary conditions (KBCs) express the 
fact that the same particles always form the free- water surface and that the fluid 
is assumed to be bounded below by a hard horizontal bed z — 0. The dynamic 
boundary condition (DBC) expresses the fact that the difference of pressure on 
the two sides of the surface 77 is balanced by the effects of surface tension. Thus, 
the boundary-value problem for capillary-gravity water waves is 



u t + uu x + vu z = -- n Vx 



v t + uv x + vv z = -~ p Pz - <l 



(EE) 



u x + v z = (MC) 

v = rj t + urix on z = h + r/(x, t) (KRO 
« = onz = 
P = Po~Ti.> on z = h + i](x,t) (DBC) 

where p(x,z,t) is the pressure, g is the constant gravitational acceleration, po 
is the constant atmospheric pressure, the parameter F(> 0) is the coefficient 
of surface tension and is the mean curvature (up to a factor 1/2) of the 
surface. For a surface defined as the function r)(x,t), the mean curvature has 
the following expression 

_L 7 i xx to) 

R (1+7,2)3/2 U 

In respect of the well-posedness for the initial- value problem for (TT]) there has 
been significant recent progress; see [UJ and the references therein. 

An important category of flows is those of zero vorticity (irrotational flows), 
characterized by the additional equation 

u z - v x = (3) 
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In what follows we will consider this type of flow. The idealization of irrotational 
flow is physically relevant in the absence of non- uniform currents in the water. 



2.2 Small-amplitude approximation of the water-wave 
problem 

We search for a linear approximation of the water-wave problem ([I])-©. We 
set the constant water density p = 1. If A > is the wavelength and a > 
is the wave amplitude, we make the following change of variables (yielding the 
non-dimensionalization of the problem; see |31j ) 



x i Y Xx, z i v h$z, r\ i — y ar), 1 1— > 

A 



u i— > V ghou, v i — y ho ^ 9 / l ° v 



/gh 



(4) 



P H' po + gh a (l - z) + gh p (5) 

where, to avoid new notation, we have used the same symbols for the non- 
dimensional variables x, z, t], t, u, v, p on the right-hand side. The non- 
dimensional pressure variable measures the deviation from the hydrostatic pres- 
sure po + gh (l — z). We now apply the following scaling of the non-dimensional 
variables 

p y- > ep, (it, v) h-> e(u, v) (6) 
avoiding again the introduction of new notation. 

Taking into account ((4]) , ([5]) and ([6]) , the two-dimensional capillary-gravity waves 
on irrotational water of finite depth are described, in non-dimensional scaled 
variables, by the following boundary-value problem: 

u t + e(uu x + vu z ) = -p x 
5 2 [v t + e(uv x + vv z )] = -p z 
u x + v, = 



v = r/t + eui] x on z = 1 + erj(x, i) 

P = 1 1- (^pr) jp^Sjrjm on z = 1 + l ) 
v = on z = 



(7) 



where we have introduced the amplitude parameter e = and the shallowness 
parameter S — It is conventional to write = S 2 W e , with W e = 
a Weber number. This parameter is used to measure the size of the surface 
tension contribution. 

By letting e —¥ 0, S and W e being fixed, we obtain a linear approximation of 
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our problem, that is, 

ut + Px = 
S 2 v t + Pz = 

Ux + v z = 

u z - 5 2 v x = (8) 
v = rjt on z = 1 

p = rj- S 2 W e r] xx on z = 1 
u = on z = 

Manipulating the first four equations of system ©, we obtain that 

v zz + S 2 v xx = (9) 

Applying the method of separation of variables, we seek the solution of the 
equation © in the form 

v(x,z,t) = F(x,t)G(z,t) (10) 

Substituting (ITU|) into the equation ^ , separating the variables and taking into 
account the expressions of v on the boundaries, that is, the fifth equation and 
the last equation in ([5]), we find 

v(x, z, t) = , t-jz sm\v(k5z)ri t (11) 
smh(fco) 

where k > is a constant that might depend on time. For the component u of 
the velocity field, taking into account (fTTj) and the fourth equation of system 
((5), we obtain 

u(x, z, t) = , cos\\{k5z)r]tx + F(x, t) (12) 

k sinh(Ko) 

where J : {x 1 1) is an arbitrary function. The components u and v of the velocity 
have to fulfill also the third equation in ©, hence, in view of pT|) and (|12[) . 

■ cosh(fc(5z)?7 t2;a; + ^LJElL = — . k<> cosh(kSz)r] t (13) 



ksmh(kS) ' dx sinh(fc(5) 

The above relation must hold for all values of x € R, and < z < 1. It 
follows that the function !F(x, t) is independent of x, therefore we will denote 
this function by F{t), and 

Vtxx + k 2 r] t = (14) 
We seek periodic travelling wave solutions; thus, for the equation (fT4|) with 



k = 2vr (15) 

we choose the following solution 

7](x, t) = cos(27r(a; — ct)) (16) 
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where c represents the non-dimensional speed of propagation of the linear wave 
and is to be determined. 

In order to find the expressions of the pressure we take into account the first 
two equations in (|5]) and the expressions of the velocity field from above. Thus, 
we obtain 

2n5c 2 

p(x, z, t) = — — — cosh(27Tf)z) cos(27r(a; — ct)) + xF'(t) (17) 

sinh(27ro) 

On the free surface z = 1 the pressure (|T7)) has to fulfill the sixth equation of 
system (|8]). Hence, in view of (fl"6|) . we get 

2ttSc 2 coth(27r(5) cos(27r(a; - ct)) + xF' {t) = (1 + 4 7 r 2 f) 2 W e ) cos(2tt(x - ct)) (18) 

The above relation must hold for all values x £ R; therefore, we get 

T(t) = constant := cq (19) 

and we provide the non-dimensional speed of the linear wave 

2 tanh(27r<5) „ 2r2rlM A / 4yr 2 r\ , / 2nho\ 
C = + ^ We) = 2^o i 1 + W J tanh ("IT j (20) 

We observe thus, that the speed of propagation of the wave varies with the 
wavelength A, with the undisturbed depth ho and with the coefficient of surface 
tension T. 

Summing up, system ([5]) has the solution 
7](x, t) = cos(27r(a; — ct)) 

u(x, z, t) = sin ^^) cosh(27T(5z) cos(27r(a; - ct)) + c 

v ( x > z ' *) = ^ts) sinh(2^z) sin(27r(x - ct)) ( 21 ) 

P( x > z > *) = sinhttvr^) cosh(27T(5z) cos(2tt(:z; - ct)) 

with c given by (|2U)) . 



3 Particle trajectories 

Let (x(t), z(t)) be the path of a particle in the fluid domain, with location 
(x(0),z(0)) := (xo,zq) at time t = 0. Below small-amplitude capillary-gravity 
water waves, taking into account (|2"TT) . the motion of the particles is described 
by the following differential system 

§ = «(*. 2 > *) = iiTTO cosh(27r^) cos(2^(* - ct)) + c 
ft = *) = iiJpF*) smh(27rfc) sin(2^(x - ct)) ^ J 

The right-hand side of the differential system (|22l) is smooth and bounded; 
therefore, the unique solution of the Cauchy problem with initial data (xq, Zq) 
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is defined globally in time. 

Notice that the constant Co is the average of the horizontal fluid velocity over 
any horizontal segment of length 1; that is, 



1 



x+l 



c = - I u(s,z,t)ds, (23) 

representing therefore the strength of the underlying uniform current (see also 
[T2"]). Thus, Co = will correspond to a region of still water with no underlying 
current, Co > will characterize a favorable uniform current and Co < will 
characterize an adverse uniform current. 

To study the exact solution of system (|22j) it is more convenient to rewrite it in 
the following moving frame 

X = 2n(x - ct), Z = 2nSz (24) 

This transformation yields 

= ^ts) cosh (Z) cos(X) + 2n(c - c) 

} <¥± — 47 " 

^ dt sinh 

I) c = c 

In this case, differentiating with respect to t, system (|25p can be written into 
the following form: 

f = - ^f/,, sin(2X) 

) at 2 smh^(27r<5) \ ' /'Ofi^ 

I = i^inh(2Z) (26) 



ig^sinh(Z)sin(X) (25) 



We denote by 

A 2 := 5 (27) 

sinh 2 (27T(5) 

We observe that A 2 as a function of 8 is decreasing and lim^o A 2 = 2w 2 c 2 , 
lim^oo A 2 = 0. Thus, 

< A 2 < 2n 2 c 2 (28) 

System (|26l) integrates at 

dt )„ 2 = ^co S (2X) + Cl (29) 



(f y = A 2 cosh(2Z) + c 2 



ci, C2 being the integration constants. Because the right-hand side of the first 
equation in (|29|) has to be bigger then zero, the constant C\ has to satisfy the 
following condition: 

ci + A 2 > (30) 
For the first equation in (f2T)| we use the substitution 



1 — v 2 2v 1 

tan(X) = y , cos(2X) = — JL , sin (2X) = , = — -dy (31) 

1 + y 1 + 2/ 1 + 2/ 
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In the new variable, the first equation in (f!29")) takes the form 

(!) 2 = ^ 2 (i-y 4 ) + C i(i + 2/ 2 ) 2 ( 32 ) 

The solution of the equation (|3"2j) involves an elliptic integral of the first kind: 

± f dV = t (33) 

J ^{c 1 ~A 2 )y 4 + 2c 1 y 2 + c 1 +A 2 

The elliptic integral of the first kind from (|33|) may by reduced to the Legendre 
normal form. In order to do this we consider the substitution 

V 2 = s (34) 

Therefore, the left-hand side in (|33[) becomes 



+ i dy = ± f ds 



v ^ + ^<> 2 - " - - 12 J 2^{c 1 -A 2 )s{s + l){s + ^) 



(35) 

Further, we introduce a new variable ip. The definition of this variable depends 

on the sign of c\ — A 2 . 

If 

ci - A 2 > (36) 
then we introduce the variable ip by 

s = tan 2 tp (37) 

and we get 

(ci - A 2 ) S (s + 1) (s + *±*) = ( Cl + A 2 )^ [1 - k 2 sin 2 v ] 
\ ci — A* 1 J cos 4 ip L J 

tan tp 

as = 2 - — a<p 

cos^ ip 

where the constant < k\ < 1 is given by 

OA 2 

* = (38 » 

Therefore we obtain the Legendre normal form of the integral in (I33|) : that is, 

+ -^— f —^ =t (39) 

+ Al J x /l-fc 2 sinV 



The inverse of the integral in (|39|) is the Jacobian elliptic function sine amplitude 
(see, for example, pQ), an odd periodic function of order two, 

sn (±y/ci + A 2 t:k^j :=sin<^ (40) 
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In view of the notations ([3"T]) . we get that 

, , sn ( \l c\ + A 2 1: k\] / , \ 

y(t) = ± yJ—^ '( := ±sc (V^T^i; h) (41) 

cn (yci + A 2 t; fci J V / 

where cn(Vci + A 2 t; fci) := cos(p is the Jacobian elliptic function cosine am- 
plitude, an even periodic function of order two, and sc is Glaisher's notation for 
the quotient sn/cn (see, for example, [T]). 
If 

ci - A 2 < (42) 

then, taking also into account the condition ()30[) . we introduce the variable </? 
by (see [38] Ch. VI, §4, page 602) 

A2 + C l 2 /.ON 

S = A 2 - ci C ° S ^ ( ^ 

In this case we get (see [30] ) 



, v4 2 - Cl 

where the constant < k 2 < 1 is given by 



y(;)=±j4J ±£i cn (V2A 2 i;fc 2 ) (4-1) 



tanh(Z) = w , cosh(2Z) = , dZ = -dw (46) 



For the second equation in ([29]) we use the substitution 

1 — w 2 1 — w 

In the new variable, the second equation in (1291) takes the form 

^ = A 2 {l-w i ) + c 2 (\-w 2 ) 2 (47) 

The solution of the equation ([47]) involves an elliptic integral of the first kind: 

/dw , „„. 

; = i (48) 

v/(c 2 - A 2 )™ 4 - 2c 2 ™ 2 + c 2 + A 2 

The elliptic integral of the first kind from (|4"8")) may by reduced to the Legendre 
normal form. In order to do this we consider the substitution 



w 2 = r (49) 



The left-hand side in (14"51) becomes 



. dw _|_ f dr 



2^{c 2 -A 2 )w±-2c 2 w 2 +c 2 + A 2 J 2 ^(c 2 -A 2 )r(r-l)(r~^§) 

"(50) 
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As in the case of the integral in (f3"5j) , we introduce a new variable tp. The 
definition of <fi depends on the sign of c 2 — A 2 and c 2 + A 2 . There are three 
possibilities: 
c 2 - A 2 > 0, 

c 2 - A 2 < and c 2 + A 2 > 0, 
c 2 + A 2 < 0. 

If 

c 2 - A 2 > (51) 
then we introduce the variable (f) by (see [38] Ch. VI, §4, page 602) 

r = sin 2 (52) 

and we get 

(c 2 - A 2 )r(r -l)(r- ^4o] = (C2 + ^4 2 ) sin 2 cos 2 (l - A; 2 sin 2 0) 



c 2 -A 2 , 
dr = 1 sin cos </><i</> 

where the constant < k 2 < 1 is given by 

k 2 = ( 53 ) 
d c 2 + A 2 

Therefore we obtain the Legendre normal form of the integral in (l4"5)k that is, 
±-=J— \-=J* =t (54) 

The inverse of the integral in (|54p is the Jacobian elliptic function sn 



±y/c 2 + A 2 t;k 3 j :=sin</> (55) 
In view of the notations (|4TJ)l . (j5"2")l . we get that 



w(t) = ± sn (y C2 + A 2 t;k 3 j (56) 

If 

c 2 - A 2 < and c 2 + A 2 > (57) 
then we introduce the variable (j> by (see [38 Ch. VI, §4, page 602) 

r = cos 2 (j) (58) 

and we get 

C 2 + A 2 \ . 2 -2 , „2 i (i ,2 -2 A 



(c 2 - A 2 )r(r - 1) ( r - - - ~ ) = 2A 2 sin 2 cos 2 (l - fc 2 sir 

V C2-A A > 

dr = —2 sin </> cos 



3 Particle trajectories 



12 



where the constant < k\ < 1 is given by 



A 2 -c 2 
2A 2 



(59) 



Therefore we obtain the Legendre normal form of the integral in (l48l) ; that is 

1 



± 



1 — fcf sin 



V2A 

The inverse of the integral in (|60p is 

sn (±%/2A2t;fc 4 
In view of the notations (|49|) . (|58]). we get that 



If 



= ± cn (V2A^t; 



c 2 + A 2 < 

then we introduce the variable ^ by (see [55] Ch. VI, §4, page 602) 

2A 2 



r = 1 



C2 



yl 2 



sin 2 i 



(60) 

(61) 

(62) 
(63) 

(64) 
(65) 
(66) 



Thus, with ([31) and (|46|) in view, the solution of system ([29]) has the following 
expression 

X(t) = arctan [y(t)} 

Z{t) = arctanh [«;(«)] - ± In ±±^g {b '> 

with given by (j41]) or (j44|) and w{t) given by (J56j) or (J62J) or ((65]). From 
(f24]) and (|67| . the solutions of system (f22]) with the constant Co equal to the 
speed of propagation of the linear wave c have the following expressions: 



In this case we get (see [3"U] ) 



w(t)=±Jl- 



2A 2 



A 2 



C-2 



sn 2 (V-4 2 -c 2 ^;fc 5 ) 



where the constant < fcf < 1 is given by 

fc 2_ 2A 2 
5 ~^ 2 -c 2 



i(i) = ct ± ^arctan sc f^ci + ^4 2 t; j^rj 
z(£) = ±2^j arctanh sn ^\/c 2 + A 2 



C2+A? 



(68) 
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x(t) = ct ± T^arctan 
z(t) = ±— ^arctanh 



sc 
cn 



(Vci + A? , 



2A Z 



2A 2 



x(t) = ct± ^arctan sc (Vci + A 2 t; ^tjit ) 
z(t) = i^ijarctanh 



1 - 



A 2 +ci 



TP - 



2A 2 t . ci+Al 



x(t) = ct ± ^arctan 
z(t) = i-^arctanh sn (^TA?t; f^ 2 ") 

= ct ± ^arctan 



2T4 2 "" 



A 2 +c 
A 2 -c 



z(t) = i^^arctanh 



/o 42 +. ^4 -eg 
z ^ 2A 2 



(69) 



(70) 



(71) 



(72) 



a;(t) = ct ± ^arctan 
z(i) = i^ijarctanh 



A 2 + 

A 2 — 



^ cn(\/lL4t;^ 



2A 2 



2A 2 
A 2 -c 2 



2A 2 

' A 2 -C 2 



(73) 



Let us analyze in more detail the solution (1681) . Taking into account the ex- 
pressions for the derivatives of sine amplitude and cosine amplitude (see, for 
example, [1]), that is, 



dt 
d_ 

dt 



sn (t; k) = cn (t; fc) dn (t; fc) 
cn (i; fc) = — sn (t; fc) dn (t; fc), 



where 



dn (t; fc) := y/l - k 2 sn 2 (t;k), 
we get the derivatives with respect to t of x(t) and z(i) from (j6"5)l : 



= c ± ^Edn (V^T^t; ^r) 



2^ cn(y^t;^) 



(74) 
(75) 

(76) 



For the alternative with "+" in the expression (|76l) of x'(t), we obtain that 
x'(t) > 0, for all t. The sign of the derivative z'(t) from (f76|) depends on the 
sign of the periodic Jacobian elliptic function cn. Thus, we get, for example, 



x'(t) > 0, z'(t) < for cn Uc 2 + A 2 t; f=C ) < 



x'(t) > 0, z'(t) > for cn ( \J c 2 + A 2 t; 



c 2 +A 2 

c 2 -A 2 
c 2 +A 2 



> 



(77) 
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In this case, is the particle trajectory an undulating curve to the right? 
We observe that for that i's, denoted t + K, with K a period, for which the 
periodic Jacobian elliptic function sn (^Jc 2 + A 2 t; ) = ±1, we have 

lim :c(t) = finite := x + K , lim z(i) = ±oo (78) 
t->(t+i<") s-s-(t+iT) 

Therefore, at x = x + K the graph of the curve ([68jl will be asymptotic. 
Using Mathematica, for example, for 5 = 1, c = 10, by ([27J) we get A 2 = 1.08704, 
and choosing a = 7.91296 > A 2 , c 2 = 2.91296 > A 2 , the graph of the curve 
Q68p with "+" in the expressions of x(t) and z(i), is drawn in Figure 1. For 
5 = |, c = 10, from (gfj we get yl 2 = 146.07, and with ci = 177.93 > A 2 , 
c 2 = 253.93 > A 2 , the graph of the curve (|68p with "+" in the expressions of 
and z(t), looks like in Figure 2. 



Figure 1 




Figure2 




For the alternative with "-" in the expression (f?6"|) of a;'(t) we obtain that 

if ci + yl 2 < 4tt 2 c 2 then x'{t) > for all t 
if ci - A 2 > 4tt 2 c 2 then a;'(t) < for all f 
if a - A 2 < 47r 2 c 2 < ci+ A 2 then 



< 



Ci+A 2 — 47T 2 C 2 

2A 1 



x'(t) < for |sn (Vci + A 2 i; | 
x'(t) > for |sn (V^T^t; ^) | > y/ ^^ 

(79) 

The sign of the derivative z'[t) from (J7SJ) depends on the sign of the periodic 
Jacobian elliptic function cn; that is, the sign of the derivative z'(t) alternates 
successively, z'(t) < and z'(t) > 0. 

In this case, is the particle trajectory (|6"B"|) an undulating curve to the right, or 
an undulating curve to the left, or a looping curve? 

Using Mathematica, for 5 = ±, c = 10, A 2 = 146.07, a = 3822.93 > A 2 which 
satisfies a - A 2 < Air 2 c 2 < c\ + A 2 and c 2 = 2353.93 > A 2 , the graph of the 
curve (|6"5|) with "-" in the expression of x(t) and "+" in the expression of z(t), 
is drawn in Figure 3. 

Figure3 




Using Mathematica, one can also draw the other solutions (f69j) -(f73 |) . For 
example, for S = ±, c = 10, A 2 = 146.07, -A 2 < a = 46.07 < A 2 and 
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C2 = —253.93 < —A 2 , the graph of the curve ([73")) with "+" in the expressions 
of x(t) and z(t), is drawn in the figure below. 

Figure4 




We remark that the curves obtained are not closed curves. 
II) c ^ c 

Differentiating system (|25|) with respect to t, we get 

— + &tan(X) — h A 2 sin(2X) - b 2 tan(X) = 

at 1 at 

where A 2 is the constant from (|2T1) and 

b := 2tt(cq — c) 

Using the substitution (|31l) . the equation (ISUl) takes the form 



rff 2 



i + y 2 v * 



dt 



(80) 



(81) 



(82) 



This differential equation can be written as an Abel differential equation of the 
second kind (see [3U] )• It is solvable and its solution has the parametric form 
(for more details, see [50]): 



y(r) = ±, 



2A 2 



1, 



(C-61n|r + Vr 2 — 2^4 2 |) 
C is a constant, and the relation between t and r is the following: 

1 

f 



: dr 



VT 2 -2A 2 yV 2 - 2A 2 - (C - bin |r + Vr 2 - 2A 2 |) 2 
Thus, taking into account (|31[) . we obtain 

X(t) = arctan [y(t)} 7 



(83) 



(84) 



(85) 



with y(r) given by (1831) and r given implicitly by (|84|) . 

In order to determine Z{t) from system (|25p . with (|85p in view, we write the 
second equation of this system in the form 

dZ 4n 2 Sc . . 47T 2 r)c y(t) , ,„„. 

= ; sin(arctan w(t)l) = ; , = = dt (86) 

sinh(Z) sinh(27r5) V L!A n ' sinh(27r<5) 



If 



4ir 2 5c 



sinh(27r<5) y/l + y*{t) 



dt + const < 



(87) 
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then we get (see [50] ) 

Z(t) = 2arctanh 



exp 



Att 2 5c 



y ^ dt + const | 

sinh(2^5) y/l + y*(t) J 



(88) 



From (j2"4")l , (|85l) and (1551) . the solution of system (|22l) is written now as 
x(t) — ct + T^arctan [y(t)] 



z(t) = ^arctanh 



exp ( f , 4 .^ »W dt + const ] 



(89) 



with y(r) given by (|83p and r given implicitly by (|84[) . Taking into account 
and (J57J), we get 



47T 2 (5c 



h(27T(5) y/1 + y*(t) 



dt 



1 r(t) - V2A 
- lor 



sin 



r(t) + V2A 



(90) 



Thus, the solution in (|8"9"|) has, as function of the parameter r, the following 
form: 



x(t) = cj 



VT' ! -2A 2 y / T 2 -2A 2 ~(C-bln\T+VT 2 -2A 2 \y 

±77-arctan 



dr 



r 2 ~2A 2 



(C-bln\T+Vr 2 -2A 2 \y 



- 1 



(91) 



z(t) = arctanh 
We also remark that the above curve is not a closed curve. 



T-V2A 
T+V2A 
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